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Pairing symmetry is important to indentify the pairing mechanism. The analysis becomes par- 
ticularly timely and important for the newly discovered iron-based multi-orbital superconductors. 
From group theory point of view we classified all pairing matrices (in the orbital space) that carry 
irreducible representations of the system. The quasiparticle gap falls into three categories: full, 
nodal and gapless. The nodal-gap states show conventional Volovik effect even for on-site pairing. 
The gapless states are odd in orbital space, have a negative superfluid density and are therefore 
unstable. In connection to experiments we proposed possible pairing states and implications for the 
pairing mechanism. 



Introduction The newly discovered family of iron- 
based ReOFeAs(Re = La, Ce, Pr, etc.) high temper- 
ature superconductors are raising great interests in the 
community. [l| The superconductor consists of layers of 
FeAs which is believed to be the conducting planes. 
The ReO layers in between the FeAs layers stabilize 
the structure and donate carriers to the FeAs layers. 
For example, substitution of O by F or introducing O 
vacancies [2] dopes electrons into the system, and con- 
trarily substitution of Re by alkaline-earth elements can 
realize hole doping. [3] Some preliminary experimental re- 
sults on the superconducting (SC) state have been ob- 
tained. The specific capacity measurement [4] and nu- 
clear magnetic resonance (NMR) [5[ indicate a nodal gap, 
and the existence of Andreev bound state implies a non- 
trivial phase structure of SC gap@, while both s-wave 
and dirty d-wave behaviors are suggested by penetrate 
length measurement • On the theoretical side, local- 
density-approximation (LDA) calculations show that the 
states near the Fermi level are largely contributed by the 
five d-orbitals. [8[ Moreover, some authors demonstrated 
that it is sufficient to keep only a few of them, for exam- 
ple the d xz and d yz orbit als, to reproduce qualitatively 
the LDA Fermi surface topology. [9[ As for the driving 
force for superconducting pairing;, both electronic [l^ and 
phonon- mediated mechanisms [11] are proposed, and var- 
ious pairing symmetries are predicted. 

Since the pairing symmetry is related to the pairing 
mechanism, a classification of all possible pairing sym- 
metries [HI is important. This is more so given the fact 
that the pairing function becomes an orbital- wise matrix 
function in the multi-orbital case, which we elaborate in 
this paper. The main results are as follows. 1) From a 
two-orbital (d xz + d yz ) model we classify all possible on- 
site and bond-wise pairing basis matrices. In addition 
to the momentum dependence, the matrices themselves 
carry nontrivial symmetries, so that even on-site pairing 
may lead to a nodal or gapless pairing. 2) In the gap- 
less case the density of states (DOS) at the fermi level is 
enhanced by the SC order. 3) Only the nodal-gap cases 
show Volovik effect in an applied magnetic field. 4) Most 
surprisingly, the gapless SC state may have a negative su- 



perfluid density, and is therefore unstable against phase 
twisting. 5) In connection to available experiments we 
propose possible pairing bases and further experiments 
to reduce the candidate list. 

Symmetry Analysis In the FeAs layer there are two Fe 
ions per unit cell because of the As ions. In order to sim- 
plify the analysis, we adopt a two-band model, i.e. we 
keep the two atomically degenerate d xz and d yz orbitals 
of Fe ion, which are important for superconductivity, and 
neglect d xy , d x i_ y i and d 3z 2_ r 2 orbitals of Fe ions and all 
orbitals of As ions for a moment. Thus we will focus on 
the Fe-lattice, for which there is only one Fe ion in each 
unit cell. The effect of As ions can be partially included 
in the effective hopping integrals for the d-orbitals. We 
define x- and ^/-direction unitary vectors as connecting 
the nearest neighbor Fe-atoms. The space group of our 
model is PA/mmm, which has higher symmetry than 
the space group of whole system PA/nmm. The normal 
state of the model can be described by, in the momen- 
tum space, H = ^2 ka (/){ a ^ka where <j>\ = {d\ z , d\ z ) k(T , 
and £fc = e k To + 5fc7"3 + 7/en in which To, 1,3 are unit and 
Pauli matrices defined in the orbital space. For point 
group operations under concern, the d xz and d yz orbital 
wave functions transform as x and y, respectively. This 
dictates that <p\.ro(pk transforms as Ai g , c^Tifa as B<i g <> 
<t>\T2<j>k (which is actually absent in Hq but is included 
here for later use) as and <p\rs(pk as B\ g . Equiv- 
alently we claim that the r-matrices carries the above- 
mentioned irreducible representations, without referring 
to the orbital wave functions further. Finally to leave Ho 
invariant, e^, 5 k and jk must transform as Ai g , B\ g and 
B2 g , respectively. The concrete form of the dispersions 
in Ho (see below) is irrelevant at this stage. 

We now discuss the pairing symmetry. The system 
is invariant under spin-SU(2) transformation and one 
can classify the pairing states into spin-singlet and spin- 
triplet cases, which we assign in the last step accord- 
ing to the global antisymmetry of the pairing function 
with respect to the combined exchange of spin, orbital 
and spatial position. We therefore concentrate first on 
the symmetry of pairing as a function of momentum 
and a matrix in the orbital space. Since the pairing 
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FIG. 1: (Color) Minimal quasiparticle excitation energy as a 
function of momentm in the normal state (a), a full-gap SC 
state (b), a nodal-gap SC state (c) and a gapless SC state (d). 
Thick white lines (symbols) indicate the Fermi surface (nodal 
points). 



operator <j> aa (-k)T?P<l)p a '(k) (where <j> a ^ = d XZiVZ and 
i = 0, 1, 2, 3) transforms identically as (/>a a (k)T^ (/)/3 a f (k) , 
we immediately see that the r-matrices in the pairing ma- 
trix transform exactly as they do in Ho. For on-site pair- 
ing, ro, 1,2, 3 carries the irreducible representations. Since 
?o, 1,3 are even in orbital space and transform as Ai g , B 2g 
and Big under the point group, the spin channel must 
be a singlet. On the other hand, ir 2 is odd in orbital 
space and transform as A 2g , the spin channel must be in 
a triplet. 

The extension to pairing on bonds is almost straight- 
forward. We may multiply the above mentioned r- 
matrices by trigonometric basis functions to form the 
pairing matrices corresponding to pairing on bonds in 
real space. For the system under concern, the basis func- 
tions for nearest-neighbor bonds are cos k x + cos k y (Ai g ), 
cosk x — cos k y (Big), (sin k x , sin k y ) (carrying the two- 
dimensional Eg representation). The symmetry of the 
product is easily seen from that fact that the r-matrices 
only carry one-dimensional representations. For exam- 
ple, Big -Ti ~ Big • B 2g ~ A 2g . Since the orbital parity is 
even, the spin channel must be a singlet. In Tab. I we list 
all possible pairing basis matrices for on-site pairing and 
nearest-neighbor-bond pairing (extension to longer-range 
pairing is trivial), together with the irreducible represen- 
tations they carry, the spin symmetry, orbital parity, and 
the behavior of the corresponding quasiparticle excita- 
tion gap in the momentum space (which will be discussed 
later). We notice that even on-site pairing can carry a 
non-trivial representation, indicating the unique role of 
the two (atomically) degenerate d-wave orbitals. 



TABLE I: Pairing basis matrices carrying irreducible repre- 
sentations of the model. The first column is the index num- 
ber, the second and the third columns list the representations 
and the basis matrix functions. The spin and orbital pari- 
ties are shown in the forth and the fifth columns. The last 
column describe the behavior of quasiparticle excitation gap 
in the momentum space. Notice that To, 1,2,3 transform as 
Ai g , B2g,A2 9 , Big respectively due to two xz and yz orbitals. 
These rules are important to identify the global symmetry 
and the gap behavior of the basis matrix function. 



No. 


IR 


Basis 




Spin Orbital Parity 


Gap 


1 


Alg 


TO 




S 


E 


Full 


2 


Alg 


(cos k x + cos ky] 


TO 


S 


E 


Nodal 


3 


Alg 


(cos k x — cos k y ] 


T3 


s 


E 


Nodal 


4 


A2g 


(cos k x — cos k y ] 


Tl 


s 


E 


Nodal 


5 


Big 


T3 




s 


E 


Nodal 


6 


Big 


(cos k x — cos ky] 


TO 


s 


E 


Nodal 


7 


Big 


(cos k x + cos k y ] 


T3 


s 


E 


Nodal 


8 


B 2 g 


Tl 




s 


E 


Nodal 


9 


B 2 g 


(cos k x + cos ky] 


Tl 


s 


E 


Nodal 


10 


Eg 


J sink x iT2 
I sin k y iT2 




s 


O 


Gapless 


10' 


Eg < 


\ (sin k x + i sin k y 
[ (sin k x — i sin k y 


)ir 2 
)ir 2 


s 


O 


Gapless 


11 


A2g 


IT 2 




T 


O 


Gapless 


12 


A2g 


(cos k x + COS k y )lT2 


T 


O 


Gapless 


13 


B 2 g 


(coskx — cos ky) 


IT 2 


T 


O 


Gapless 


14 


Eg 


( sinkxTo 




T 


E 


Nodal 



14' E c 



15 E Q 



15' E c 



16 E a 



16' E, 



{ 
i 
i 



(sin k x + i 
(sin k: 



sin k y To 

sin k y )ro 

! T 



— ism 
sin k x r 3 
sin k y T3 
kx ~~\~ 1 sin ky 
kx i sin ky 
sin k x Ti 
sin k y n 

sin k y )n 
sin k y )ri 



(sin 
sin 



'sin k x + i 
sin k : 



'T 3 
'T 3 



T 
T 
T 
T 
T 



E 
E 
E 
E 
E 



Full 

Nodal 
Nodal 
Nodal 
Nodal 



Quasiparticle excitation gap To illustrate the concrete 
gap structure in the various SC states, we need the 
band structure of iron-based superconductors. Here we 
adopt the tight-binding model introduced by Ref.[13]. 
The BdG Hamiltonian is given by H = ^2 k ^\H^k^ 
where we recall that the momentum is defined in the 
"large" Brillouin (BZ) corresponding to one Fe per unit 

cell, ijj k (d X zMi d yz,kh d lz-ki> d lz-ki) T is the four_ 
component spinor in the orbital and Nambu space, and 
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Hh is, in a form of block-matrix, 
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Here = e k TQ -\-SkTs + 7/cTi is a tight-binding dispersion 
defined in i^o, with e& = — (ti + t2)(cosk x + cosk y ) — 
4ts cos k x cos k y — fi, 5k = — (ti — t2)(cosk x — cos/c^), 
and 7fe = —4^4 sin fc^ sin fc^. Here ti = —1, £2 = —1.3, 
t 3 = t 4 = -0.85, and /i = 1.45, in unit of \ti\. The 
gap amplitude V = 0.4 is chosen for illustration, and 
Afc is selected from the basis matrix functions listed 
in Tab. I. The hamiltonian can be exactly diagonalized 
to obtain the quasiparticle excitations and other SC 
properties. Zero energy excitations are determined by 
det(Hk) = 0. In the case of SC state No.l in Tab. I, 
det(H k ) = (e 2 -S 2 -j 2 - V 2 ) 2 + (2Ve k ) 2 = has no so- 
lution. This corresponds to the full-gapped case. For SC 
state No.5, det(iTfc) = (e 2 - 5 2 - 7 | + V 2 ) 2 + (2VS k ) 2 = 
is satisfied at four sets of nodal points in BZ. Another in- 
teresting example is SC state No.ll, for which det(Hk) = 

(4 - s l - il + y2 ) 2 = holds alon g lines in BZ - Such 

zero-lines form the fermi surface (FS) (slightly different 
from the normal state FS) for the BdG quasiparticles. 
This is the gapless case. 

Due to the two-orbital character of the SC state, the 
nodal points or the BdG FS may be located away from 
the normal state FS. We therefore need the quasiparticle 
excitation energy in the whole BZ to characterize the 
gap structure. Several typical results are shown in FigO] 
where the minimal quasiparticle excitation energy as a 
function of momentum in (a) the normal state, (b) a 
full-gap SC state No.l, (c) a nodal-gap SC state No.5 
and (d) a gapless SC state No.ll. Here thick white lines 
(symbols) highlight the zero excitation energy contour 
(nodal points). By checking the excitation spectra of all 
possible cases in Tab. I, we summarize that: (i) All t\^- 
bases are nodal, consistent with the fact that these r 
matrices carry £>2p,i# representations, (ii) All T2-bases 
(which carries odd orbital-parity) are gapless (unless the 
pairing energy scale is of the order of the band width), 
(hi) The ro-bases leads to full or nodal gaps, depending 
on the momentum basis function. 

Volovik effect and superfluid density In this section, we 
discuss the physical properties of various SC states, espe- 
cially the DOS, Volovik effect and superfluid density. We 
did the calculations for all SC states listed in Tab. I. Since 
we found that the qualitative behavior is the same for SC 
states belong to the same category (full-gap, nodal-gap 
or gapless) , we only present the results for representative 
states, namely, the full-gap SC state No.l, the nodal- 
gap SC state No.5 and gapless SC state No.ll. Fig|2|(a) 
shows the DOS. The low energy U-shaped (V-shaped) 
DOS in the full-gap (nodal-gap) SC states are conven- 
tional. The gapless case is rather exotic. The low energy 
DOS is not gapped but slightly piled up by SC order, 
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FIG. 2: (Color) Physical properties of various SC states: 
(a) Density of states; (b) Change of zero energy DOS by a 
magnetic field H oc q 2 ; and (c) Low temperature superfluid 
density. 



due to the presence of BdG FS. FigEJb) shows the effect 
of a magnetic field on the DOS, which is calculated by 
a semi-classical method as follows, [li) The effect of vor- 
tices is simulated by averaging over the direction of the 
superflow momentum q with q oc y/~H. The change in 
zero energy DOS A7V(0, q) = (JV(0, q)) - 7V(0, 0) can be 
probed by the magnetic field dependent specific heat. [4] 
It is seen in Figj2jb) that the Volovik effect is absent 
in the full-gap SC state, but is manifest as a linear rise 
of AN with q in the nodal-gap SC state. Moreover, for 
the gapless state, the presence of magnetic field actu- 
ally reduces the DOS (as seen for moderate q). FigJ2fc) 
shows the temperature dependence of superfluid density. 
In the low temperature limit we neglect the temperature 
dependence of the pairing amplitude. The full-gap and 
nodal-gap SC states exhibit activated and linear drop, 
respectively. Contrarily, the superfluid density is nega- 
tive at all temperatures in the gapless case. This arises 
from the fact that piling up of low energy DOS leads to 
an overwhelming paramagnetic contribution against the 
diamagnetic part. It indicates that the gapless state is 
unstable upon phase twisting, implying a tendency to- 
ward possible Fulde-Ferrell-Larkin-Ovchinnikov state or 
magnetic ordering. 

Discussion Basing upon our analysis presented above 
and recent experiments, we could discuss the possible 
pairing symmetry of iron-based superconductors. Spe- 
cific heat [4], NMR[5] penetrate length [7] and tunneling 
[6[ measurements are consistent with spin-singlet nodal 
pairing with sign changes. On the other hand, on-site 
singlet pairing is unfavorable to Hund's rule coupling, 
while bond-wise triplet pairing is inconsistent with the 
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antiferromagnetic exchange. [15j We are then left with the 
bond-wise singlet pairing cases No. 3, 4, 6, 7, and 9 in Tab. I. 
Out of these cases we observe that only the A2 g case 
No. 4 and the I>2g case No. 9 have nodal points in the x- 
and ^-directions, which are more relevant if the electron 
pockets are important. We propose that Raman scatter- 
ing and phase-sensitive probes could further reduce the 
redundancy. 

Some remarks are in order before closing. First, it 
should be emphasized that the symmetry classification 
is robust, while the assignment of a particular pairing 
symmetry has to be sharpened or even altered accord- 
ing to further systematic and intrinsic experimental re- 
sults. Second, a general pairing matrix can be decom- 
posed into a linear combination of the bases. In prin- 
ciple, bases belonging to different irreducible represen- 
tations do not mix, but mixing in other cases can not 
be ruled out on symmetry ground alone. In particular, 
the two-dimensional representations, like the case No. 10 
may be mixed into one of the chiral-symmetry-breaking 
states shown in No. 10'. Third, even if the effect of As 
ions neglected so far are reconsidered, the point group is 
still D^, so that the above symmetry classification still 
holds. Finally, several authors propose models including 
d xy or even all five d-orbitals.[16j These may be included 
in our analysis and result in more pairing bases. For ex- 
ample xy can pair with yz, which transform as x, forming 
one component of the E representations. 

We became aware of a related work [17] after we fin- 
ished the present paper. 
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